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ABSTRACT 


High  quality  long  fiber  reinforced  composites,  such  as  those  used  in 
aerospace  and  industrial  applications  are  commonly  processed  in  autoclaves. 

An  adequate  resin  flow  model  for  the  entire  system  (laminate/bleeder/ 
breather),  which  provides  a  description  of  the  time-dependent  laminate 
consolidation  process,  is  useful  in  predicting  the  loss  of  resin,  heat 
transfer  characteristics,  fiber  volwne  fraction  and  part  dimension,  etc., 
under  a  specified  set  of  processing  conditions.  This  could  be  accomplished  by 
properly  analyzing  the  flow  patterns  and  pressure  profiles  inside  the  laminate 
during  processing. 

In  this  paper  a  newly  formulated  resin  flow  model  for  composite  prepreg 
lamination  process  is  reported.  This  model  considers  viscous  resin  flows  in 
both  directions  perpendicular  and  parallel  to  the  composite  plane.  In  the 
horizontal  direction,  a  squeezing  flow  between  two  nonporous  parallel  plates 
is  analyzed,  while  in  the  vertical  direction,  a  poiseuille  type  pressure  flow 
through  porous  media  is  assumed.  Proper  force  and  mass  balances  have  been 
made  and  solved  for  the  whole  system.  The  effects  of  fiber-fiber  interactions 
during  lamination  are  included  as  well.  The  unique  features  of  this  analysis 
are  (i)  the  pressure  gradient  inside  the  laminate  is  assumed  to  he  generated 
from  squeezing  action  between  two  adjacent  approaching  fiber  layers,  and  (ii) 
the  behavior  of  fiber  bundles  is  simulated  by  a  Finitely  Extendable  Nonlinear 
Elastic  (FENE)  spring.  Favorable  comparisons  between  model  predictions  and 
experimental  data  available  in  literature  are  found. 


f**®^**®  page  blank  not  nLMtu 


FOREWORD 


This  work,  performed  at  NASA  -  Langley  Research  Center  (LaRC),  was 
Supported  by  the  LaRC  Polymeric  Materialr  Branch  (PMB)  under  Contract 
NASl-18000.  This  is  a  progress  report  on  an  ongoing  research  project  toward 
monitoring  and  controls  of  composite  laminate  fabrication  processing  in 
autoclave.  R.  M,  Baucom  (PMB)  was  the  Technical  Monitor.  This  work  also 
represents  an  extended  version  of  the  work  entitled  “A  Theoretical  Study  of 
Resin  Flows  for  Thermosetting  Materials  During  Prepreg  Processing"  previously 
reported  by  NASA  CR-172442,  July,  1984.  Thanks  are  due  to  James  Shen  (ODU) 
who  helped  to  check  numerical  solutions  of  Eqs.  (9)  and  (11)  as  plotted  in 
•  igure  5;  to  Dr.  .1,  A.  Hinkley  (PMB)  who  reviewed  the  entire  manuscript  prior 
to  1t^  publication. 
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AUTOCLAVE  PROCESSING  FOR  COMPOSITE  MATERIAL  FABRICATION 

I.  AH  ANALYSIS  OF  RESIN  FLWS  AND  FIBER  COMPACTIONS  FOR  THIN  LAMINATES 

I.  INTRODUCTION 

High  quality  long  fiber  reinforced  composites,  such  as  those  used  in 
aerospace  and  industrial  applications,  are  commonly  processed  in  autoclaves. 
During  processing,  the  composite  materials  are  subjected  to  prescribed 
elevated  temperatures  and  pressures.  Selection  of  a  cure  cycle  (i.e., 
temperature  profile)  will  dictate  the  kinetics  of  the  reactive  resin  matrix, 
and  consequently  the  profile  of  chemoviscosity  built-up.  Increased  molding 
pressure  during  processing  will  result  in  resin  flows  both  perpendicular  and 
parallel  to  the  fibers.  The  applied  pressure  helps  to  consolidate  the 
composite  laminates  and  to  squeeze  out  excess  resin  and  voids.  An  adequate 
resin  flow  model  for  the  entire  system  (laminate/bleeder/breatner),  which 
provides  a  description  of  the  time-dependent  laminate  consolidation  process, 
is  useful  in  predicting  the  loss  of  resin,  heat  transfer  characteristics, 
fiber  volume  fraction  and  part  dimension,  etc.,  under  a  specified  set  of 
processing  conditions.  This  could  be  accomplished  by  properly  analyzing  the 
flow  patterns  and  pressure  profiles  inside  the  laminate  during  processing. 

Considerable  work  has  been  conducted  in  the  past  by  numerous  researchers 
in  searching  for  a  relationship  between  chemoviscosity  and  cure  kinetics 
[1,2].  However,  the  mechanics  governing  the  flow  of  resin  associated  with  the 
composite  lamination  process  has  received  little  attention.  Springer  and  Loos 
[3-6]  and  Lindt  [7]  had  presented  pure  viscous  flow  models.  In  their 
analyses,  the  applied  pressure  during  composite  prepreg  lamination  process  is 
carried  by  the  resin  matrix  only.  The  pressure  gradients  resulted  inside  the 
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laminate  create  flows  In  both  directions  perpendicular  and  parallel  to  the 
planes  of  the  composite.  Although  favorable  comparisons  between  model 
predictions  and  experimental  observations  have  been  reported  for  some  cases, 
it  is  generally  recognized  that  an  adequate  flow  model  for  the  composite 
consolidation  process  must  include  the  effect  of  compacted  fibers  as  well. 
Recent  measurements  [8]  indicated  that  resin  pressure  could  be  as  low  as  1/3 
to  1/10  of  the  applied  pressure,  and  the  missing  load  must  be  carried  by  the 
fibers.  Bartlett  [9]  developed  a  theoretical  flow  model  for  glass-reinforced 
resin  during  lamination  of  multilayer  printed  circuit  boards  in  the  electronic 
industry.  His  analysis  appeared  to  be  the  first  attempt  to  take  into  account 
the  effects  of  layers  of  glass  fabric  coming  into  contact  with  one  another  as 
the  resin  is  squeezed  out.  The  model  was  later  compared  with  experimental 
results  by  Bloechle  [10,11]  using  a  paral lel-plate  plastometer,  and  was 
successfully  applied  to  quality  control  of  Incoming  epoxy  B-stage  prepregs  in 
a  manufacturing  environment. 

More  recently  Gutowski  [12]  reported  a  resin  flow/^iber  deformation  model 
for  composites.  By  assuming  fibers  possess  small  curvatures,  an  elastic  fiber 
model  depicted  by  a  rapidly  stiffening  spring  was  constructed.  The  model  was 
also  shown  to  fit  experimental  data  favorably. 

In  this  paper  a  newly  formulated  resin  flow  model  for  composite  prepreg 
lamination  process  is  reported.  This  model  considers  viscous  resin  flows  in 
both  directions  perpendicular  and  parallel  to  the  composite  plane.  The 
effects  of  fiber-fiber  interactions  during  lamination  are  included  as  well. 

The  unique  features  of  this  analysis  are  (i)  the  pressure  gradient  inside  the 
laminate  is  assumed  to  be  generated  from  squeezing  action  between  two  adjacent 


approaching  fiber  layers,  and  (ii)  the  behavior  of  fiber  bundles  is  simulated 
by  a  Finitely  Extendable  Nonlinear  Elastic  (FENE)  spring.  Comparisons  between 
model  predictions  and  experimental  results  will  also  be  made. 
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2  dp/dz  (psf/inch) 

Viscous  flow  Activation  Energy  (Kcal/mole) 

Viscosity  cure  Activation  Energy  (Kcal/mole) 

External  load  (lbs) 

Separation  (inch)  between  parallel  plates 
Initial  separation  (inch)  between  parallel  plates 
2  dh/dt  (inch/sec) 

Permeability  (inch2)  of  porous  material 
Viscosity  rate  constant 
Material  constant  (min-*) 

Ambient  pressure  (psi) 

Pressure  (psi)  absorbed  by  glass  fabric  or  fiber  bundles 
Defined  by  Eqs.  (12),  (13) 

Universal  gas  constant  (kcal/mole  ®K) 

Cylindrical  coordinate  system 
Curing  temperature  (®K) 

Curing  time  (sec) 

?n’p5JoJs^la"feMlf’  of  space 

Velocity  (inch/sec)  in  r  and  z  direction,  respectively 
^esin'nJwr^^"  tnickness  (inch)  of  porous  material  where 

Chcmoviscosity  (poise) 

Initial  viscosity  (poise)  at  t  «  0 
Materia’  constant  (poise) 

Oer.sity  (Ih/ft^)  averaged  over  a  region  in  porous  material 
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II.  THEORY 


Physical  Model 

A  simplified  idealized  schematic  diagram  of  multi-layer  composite 
laminates  is  shown  In  Fig,  1.  A  stack  of  five  layers  of  prepreg  tapes  Is 
confined  between  two  steel  plates  and  layers  of  porous  materials.  The  porous 
materials  Include  those  commonly  called  bleeder  and  breather  materials.  In 
reality,  the  fibers  are  not  perfectly  straight  and  aligned  as  shown.  They  are 
rather  curved  and  have  Indefinite  number  of  points  of  contact  In  between. 

When  the  laminate  compactions  occur  due  to  the  external  load,  the  number  of 
points  of  contact  among  layers  of  fiber  bundles  increases  rapidly,  and  so  does 
the  fraction  of  external  load  carried  by  the  fiber  bundles.  Such  behavior  can 
be  conveniently  simulated  by  a  non-Hookean  elastic  spring. 

In  order  to  develop  a  viscous  flow  model  for  the  resin,  a  ‘low  channel  Is 
postulated  and  shown,  schematically  in  Figure  2(a),  where  resin  is  confined 
between  two  parallel  porous  plates  separated  by  a  distance,  2h.  Upon 
application  of  a  force,  F,  the  resin  Is  squeezed  outward  horizontally  and  also 
vertically  through  the  porous  media.  It  is  assumed  that  these  two  flow 
directions  can  be  decoupled  conceptually  as  shown  in  Fig.  2(b)  and  2(c). 

Fig.  2(b)  Illustrates  a  squeezing  flow  between  two  non-porous  plates  separated 
by  a  same  distance,  2  h,  A  vertical  pressure  flow  through  porous  materials  is 
illustrated  in  Fig.  2{c)  (assuming  no  horizontal  flow),  and  that  the  flow  is 

^••iven  by  a  pressure  drop  (p-pa)  across  a  characteristic  distance  zo  where 
p  is  the  pressure  generated  by  the  squeezing  action  between  two  approaching 
piates,  and  pj  is  the  ambient  pressure.  A  mass  balance  between  change  of 

plates  separation  and  the  flows  in  these  two  directions  can  be  established. 
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The  externaHy  applied  force  F  must  be  balanced  by  the  pressures  and  the 
elastic  force  of  fiber  bundles  as  well. 

Mathematical  Formulation 

-Vertical  Flow  Through  Porous  Media 

For  toe  flo-  of  e  flufo  through  a  porous  neoium.  the  eguatlon  of  ™tloo 
can  be  replaced  by  Darcy's  law  [13] 

Vo  -  -|(FP  -  09).  ,,, 


-here  the  unoerltoed  ouaotittes  denote  sectors.  It  ,s  the  penseahillty  of  the 
porous  medlui.  „  Is  the  siscoslt,.  V,  is  a  superficial  .elodt,  averaged  over 
a  snail  region  of  space,  and  o  and  p  are  density  and  pressure,  respectl.lty, 

averaged  o,er  a  region  available  to  flo-  that  Is  large  »lth  respect  to  the 
pore  site. 


For  an  incompressible  liquid  and  constant  K  and 
the  equation  of  continuity  can  be  reduced  to 


n.  Eq.  (1)  together  with 


AS  a  first  order  approximation,  we  assume  that  Eq.  (2)  is  applicable  to  our 

system  (Fig.  2(c))  in  a  unidirectional  flow.  Eqs.  (1)  and  (2)  can  then  be 
conib^ned  to  give 


Vo 


K 

n 


'2‘ 


(3) 
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j  where  «  dp/dz  is  a  constant.  Eq.  (3)  has  also  been  adopted  by  Springer 

i  .  *  ' 

I  j  and  Loos  [3-6]  in  their  analysis. 

i 

I  .  ' 

Horizontal  Squeezing  Flow 


We  now  consider  the  squeezing  flow  between  two  non-porous  plates,  as 
shown  in  Fig.  2(b).  A  cylindrical  coordinate  system  (r.e.z)  is  chosen  for 
convenience.  The  velocity  and  pressure  fields  are  assumed  Vp  * 

■  V^Cz)  and  p  ■  p{r)  only.  The  equation  of  continuity  becomes 


h 

2  /  dz  ♦  r  (fi  ♦  Vo)  .  0,  (4, 

0 

where  fi  »  dh/dt  is  a  measure  of  the  speed  of  approach  of  the  two  parallel 
plates  upon  the  application  of  force  F,  and  »  fi  at  z  =>  h. 


For  the  velocity  field  assumed.  Bird  et  al.  [14]  have  shown  that  the 
r-ccmponent  equation  of  motion  can  be  integrated  to  give 


h^r  ,  1 

“2“ 


nr  dr 


n  -  (— 

dr'  'h'  i 


(5) 


Substituting  Eq.  (5)  into  Eq,  (4),  we  have 


i£. .  3  ^  , 

dr  7  *0'* 


(6) 
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which  can  then  be  integrated  to  obtain  the 


pressure  distribution  as 


p(r).p^  •  -  I  n  ^  (l!  *  V.)  [1  .  (J)"] 

Eq.  (7)  represents  the  squeezing  generated  pressure  drop  between  two  parallel 
plates  approaching  with  a  speed  fl(t).  The  pressure  distribution  is  shown  to 
be  parabolic,  and  is  dependent  upon  distance  of  separation  between  plates, 
h(t),  and  plates  approaching  speed  fi{t). 


Substituting  Eq.  (6)  into  Eq.  (5),  we  have 


V,(r.2)  (3, 

»hich  aescrtBes  tne  velocity  profile  of  flo«  in  the  horltootal  airectloo 
between  parallel  plates. 

C*  Elastic  Forces  from  Fiber  Bundles 

It  13  recognized  that  in  practice  the  prepreg  tapes  are  not  perfectly 
straight  and  aligned  within  the  composite  laminate.  The  fiber  bundles  are 
assumed  to  be  curved  initially  and  to  behave  as  an  elastic  spring.  As  the 
force  is  applied,  the  number  of  points  of  contact  in  between  increases 
rapidly,  and  the  elastic  spring  becomes  stiffen.  Such  behavior  can  be 

Simulated  by  a  Finitely  Extendable  Nonlinear  Elastic  (FENE)  spring  which  has  a 
force  law  of  the  following  form; 
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F,  ■  f,-  A 


(ho-l') 


1‘  -  (l 


<"  »h1ch  t,(t)  ts  the  distahc.  of  separatloo  ootweon  plates,  aod  h(0)  .  h.  Is 
the  Ihit.el  distance.  „  and  «.  are  nocOer  and  radios  of  prepre,  pi,  {fip,, 
tnindle)  respectively.  (h,H.S,)  then  represents  the  maxtaon  conpressihle 
..stance  achleveahle.  A  Is  the  s^ple  area,  t  Is  the  elastic  constant  and  n 
s  a  constant.  A  sprlnp  (fiber  bundles)  »1th  this  force  law  will  be  linear 
(Hootean)  for  s»all  compression,  but  will  set  stiffer  and  stiffen  as  the 
sprln,  (fiber  bundles)  Is  compressed;  furthermore,  the  spring  (fiber  bundles) 

cannot  be  compressed  beyond  (h,.p,).  because  Infinitely  large  expression 
force  is  then  required  according  to  Eq.  (9). 


fol lows: 


New  -e  can  equate  a  balance  In  forces  using  Eqs.  (2).  („.  e„o  „ 


F  ■  /  tp(r)  -  p^  r  Cjtj  .  Pf)  2  vrdr 


■R'tC.to  .  0.)  .  3  ^  , 


-here  F  is  the  applied  external  force  to  the  plates,  t,  is  a  characteristic 
.enpth  in  vertical  dl-ectlon  across  porous  media,  and  P,  is  me  pressure 
Absorbed  h,  the  fabric  or  fiber  bundles.  Physically  It  Is  noted  that  when 


layers  of  fibers  come  Into  contact  with  one  another  as  the  resin  is  squeezed 
outward,  they  begin  to  carry  a  portion  of  the  applied  load.  The  average 
pressure  applied  to  the  resin  is  therefore  the  difference  between  the  average 
applied  pressure  and  the  pressure  Pf  carried  by  the  fiber  bunoles. 

Eq.  (10)  can  be  rearranged  as: 


with 


^  fv  /• 
dt  -  -  *=2 


Q2 


h3] 


0i  =  (c^zq  +  P^) 


(11) 


(12) 


and 


^2  =  4  (13) 

The  time-dependent  laminate  thickness  h{t)  under  a  constant  load  F  can 
therefore  he  ca1cula^od  by  solving  Eq.  (U). 


10 


vr3c 


in.  Results  and  Discussion 


For  a  given  multi-layer  composite  laminate,  we  have 
V,(t)  ho 

7T“"ir(tT  (H) 

where  Vf{t),  h{t)  are  the  fiber  volume  fraction  and  laminate  thickness  at 
time  t  respectively.  Similarly  ho  »  h(0)  and  »  \t^{o)  are  the  initial 
laminate  thickness  and  fiber  volume  fraction  respectively.  Substituting  Eq. 
(14)  into  Eq.  (9),  we  have 

D  .  .  (VT(t)-vO)/v^(t) 

f  *  •'f  ^  (15) 

(v^-vO)/vf  ^ 

with  kf  -  k  hg/A,  and  denotes  maximum  flbe-  volume  fraction  obtainable. 
Experlmenta’  data  shown  In  Figure  3  were  taken  by  Gutowski  [12].  They  were 
obtained  by  compressing  6  inch  long  graphite  fiber  bundles  Impregnated  In 
light  mineral  oil  which  has  a  viscosity  of  0.05  Pas.  The  elastic  behavior 
shown  is  a  nonHookean  ••apidly  stiffening  spring.  By  selecting  V°  =  0.5, 

=  0.8,  k^  =  50,  and  n  =  1.475,  preiiictions  of  FENE  spring  model,  Eq.  (15), 
for  the  fiber  bundles  represented  by  the  solid  line,  are  shown  to  compare 
favorably  with  the  experimental  data  in  the  Figure.  During  the  selection  of 
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values  for  model  parameters,  it  is  noted  that  and  are  obtained  from 
experimental  observations,  while  kf  and  n  are  two  adjustable  variables. 

Bloechle  [10]  and  B'.rtlett  [9]  had  investigated  load/deformation 
characteristics  of  a  commercial  B-stage  1080  glass  fabric  with  a  thickness  of 
0.0036  inch  using  a  parallel -pi ate  plastometer.  Measurements  of  pressure 
carried  by  the  glass  cloth  as  a  function  of  thickness  h  are  reproduced  in 
Figure  4.  By  selecting  hg  «  0.0036  in.,  (hg-mRf)/hg  »  0,99,  n  »  6,  and 
kf  »  370,  the  FENE  model  Eq,  (9),  represented  by  the  solid  line,  is  shown  to 
compare  reasonably  well  with  the  data  points  within  the  experimental  range. 

The  Chernov is cos ity  profile  for  B-stage  1080  resin  under  isothermal  curing 
at  temperature  T  had  been  measured  by  Bloechle  [10],  and  can  be  represented  by 
a  dual  Arrhenius  expression  as  follows: 

n{t)  =  no  exp  [kt]  (16) 


with 


no  =  n^  exp  [aE^/RjT] 


k  »  k_  exp  [-aE|^/RiT] 


where  n  denotes  curing  viscosity;  no  denotes  zero  time  viscosity;  k  is 
viscosity  rate  constant.  Values  of  parameters  are  tabulated  in  Table  1. 
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Table  1.  values  of  parameters  for  B-stage  mo  epoxy  resin  and  glass  fabri, 

■  3.78  X  10“*^  poises 
k,  *  5.5  X  10’  min-* 

*  ?8  kcal/mole 
aE|^  ■  20  kcal/mole 
Rl  •  1.982  X  10-’  kcal/{mole  "K) 

K  «  9.0  X  10- *3  inch^ 

2o  ■  0.0004  inch 
R  »  2.93  inch 


The  flo,  data  obtained  by  plastoeeter  .ere  perfomwj  on  6  «  4.5  Inch 
«-sta9e  tpecteent  under  three  Uotherraal  cure  temperatures  of  140.  160.  and 
180"C.  Each  of  the  flow  test  specimens  .as  comprised  of  10  stacked  piles.  In 
order  to  maintain  a  p/a  ratio  of  1.16  psi/(so.  inch)  .nich  corresponds  to  a 
typical  NLB  lamination  procedure  utiliting  500  pst  on  an  18  .  24  inch  panel,  the 
pplied  force  F  acting  on  the  test  specimen  .as  maintained  at  045  lbs.  The 
average  B-stage  thickness  per  ply  n(tl  measured  for  three  isothermal  flo.  tes-s 
are  reproduced  in  Figure  5.  The  solid  lines  are  model  predictions  using  Eg. 

(II)  incorporating  the  FEi«  spring.  Eq.  (9).  for  fibers  bundles  as  discussed 
above,  values  of  parameters  used  in  calculations  are  also  included  in  Table  1. 

For  a  given  force  F.  the  model  seems  to  over-estimate  the  ccmpactlon  level 
of  the  laminate.  denotes  the  pressure  drop  for  flo.  in  vertical 

direction  across  a  porous  medium  .1th  a  characteristic  length  t,.  values  of 
Cat,  Increase  .ith  increasing  temperature.  As  the  top  and  bottom  bleeders 
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are  filled  with  the  resins  during  processing  as  a  result  of  vertical  flows,  it 
is  conceivable  that  the  flow  characteristic  would  change.  Uses  of  virlable 
values  of  CzZo  could  therefore  give  better  fits  between  model  and 
experimental  data.  It  Is  also  noted  that  the  range  of  measured  FENt  forces 
Pf.  Shown  In  Figure  4,  does  not  cover  the  experimental  values  of  h(t).  shown 
in  Figure  5,  wh.ere  a  reduction  of  only  10-20*  of  the  original  laminate 
thickness  been  accomplished.  The  final  laminate  thickness  can.  however,  be 
predicted  accurately  to  within  8*  at  all  three  curing  temperatures. 

Fractions  of  applied  load  carried  by  fiber  bundles  and  the  reactive  resin 

chemoviscosity  as  a  function  of  curing  time  are  plotted  In  Figure  6  and  7 

respectively.  It  Is  seen  that  a  lower  Initial  viscosity  can  be  obtained  for 

higher  Isothermal  curing  condition.  A  lower  viscosity  gives  rise  to  a 

faster  plate  approach  rate  fi{t)  at  the  early  stage  of  laminate  compaction,  and 
consequently  a  higher  rate  of  Increase  In  fractions  of  applied  load  carried  by 
the  fiber  bundles  as  shewn  in  Figure  6.  For  small  deformation  where  (PfA/F) 

<  1.5*.  the  force  response  of  FENE  spring  appears  to  be  Hookean.  As  the 
deformation  becomes  larger,  functions  of  PfA/F  vs.  t  plotted  In  Figure  6 
should  exhibit  rapidly  increasing  slopes  corresponding  to  a  rapidly  stiffening 
spring  as  discussed  before.  The  fact  that  opposite  characteristics  for  each 
curve  (I.e.  decreasing  slopes  as  a  function  of  time  in  Figure  6)  are  obtained 
indicates  that  the  majority  of  the  applied  load  is  carried  by  the  resin  matrix 

which  possesses  high  chemoviscosl ties  for  this  particular  case  under 
discussion. 
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The  increases  in  fiber  volume  fractions  Vf(t)/V®  during  processing  as 
Shown  in  Figure  8  are  calculated  from  Eq.  (14)  for  three  temperatures.  By 
knowing  the  initial  volume  fraction  of  resin.  Vj,  the  residual  resin  volume 
fraction  Vp  (t)  can  be  calculated  directly  from  h{t)  by  the  following 
equation: 


U  0  V  ®  ^  h(t) J  h(t)  (17) 

r  V 

The  pressure  distributions  p(r)  -  p^  for  the  resin  matrix  inside  the 
laminate  calculated  from  Eq.  (7)  are  plotted  in  Figure  9.  The  distributions 
are  parabolic  as  a  function  of  radius  r.  Lower  pressure  levels  as  seen  for 
180»C  are  consistent  with  the  results  shown  in  Fig.  6,  due  to  the  fact  that 
higher  portions  of  the  applied  load  are  carried  by  the  fiber  bundles.  The 
change  in  pressure  during  the  lamination  process  appears  to  be  small  for  all 
processing  temperatures  because  of  the  same  reasons  discussed  before. 
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IV,  CONCLUSIONS 


A  resin  flow  model  for  the  long  fiter  reinforced  composite  prepreg 
leminetion  process  has  teen  formulated.  The  model  considers  .iscoos  resin 
flows  in  both  directions  perpendicolar  and  parallel  to  the  composite  plane 
Tbe  effects  of  fiber.fiher  interactions  are  included  as  well.  Proper  force 
end  mass  balances  have  been  made  and  solved  foe  the  entire  system  which  is 
composed  Of  laminates,  bleeders,  and  breathers.  Loss  of  resin,  fiber  volume 
fraction,  pressure  distribution  Inside  tfe  laminate,  part  dimensions  and 
temperature  effects,  etc.,  as  a  function  of  processing  time  can  all  be 
Simulated.  By  assuming  proper  velocity  fields  for  the  sgueeting  flows  between 
parallel  plates,  a  parabolic  pressure  distribution  within  the  resin  is 
Obtained.  Tne  effects  of  fiber-flber  interactions  during  laminate  compaction 
are  simulated  by  a  Finitely  Eatendable  nonlinear  Elastic  (FENE)  force,  which 
behaves  as  a  rapidly  stiffening  spring,  and  correctly  depicts  the 
experimentally  observed  behaviors  of  fiber  responses  during  the  lamination 
process.  Comparisons  of  model  predictions  and  one  set  of  experimental  data 
found  in  the  literature  show  that  the  final  laminate  thickness  can  be 
predicted  accurately  to  within  under  various  isothermal  curing  conditions. 
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